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Abstract 

We study discrete rogue waves in an array of nonlinear waveguides. We show that very small degree of disorder due to 
experimental imperfection has a deep effect on the formation of discrete rogue waves. We predict long-living discrete 
rogue wave solution of the discrete nonlinear Schrodinger equation. 
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1. Introduction 

Rogue waves, sometimes known as freak waves 
or extreme waves, are waves that appear on a fi¬ 
nite background as a result of modulational instabil¬ 
ity. The height of rogue waves is defined as at least 
two times higher than the average surrounding back¬ 
ground. Rogue waves were observed long time ago 
in oceans. The well-known one-dimensional nonlin¬ 
ear Schrodinger equation (NTS) with attractive nonlin¬ 
ear interaction is a model equation to investigate rogue 
waves theoretically. In 1983, Peregrine found an ana¬ 
lytic solution of the nonlinear Schrodinger equation [Tl. 
The Peregrine soliton, limiting case of Kuznetsov and 
Ma soliton ||2l[3l and Akhmediev breather ||4l, explains 
how rogue waves appear from nowhere and disappear 
without a trace. The Peregrine soliton is formed from 
slightly modulated uniform background and grows until 
it reaches its maximum value at a specific time. Then 
the amplitude of the soliton decreases while the width 
increases and finally it vanishes. That is why the Pere¬ 
grine soliton is known as doubly localized wave (local¬ 
ized both in space and time). It is well known that some 
physical systems such as optics, plasma and ultracold 
atoms are also described by the nonlinear Schrodinger 
equation. Therefore, the existence of the Peregrine soli¬ 
ton, or more generally rogue wave, is not restricted to 
oceans 0161121. The experimental realization of Pere¬ 
grine soliton was first made in an optical system in 2010 
m and then in a water wave tank in 2011 g). These ex¬ 
perimental realizations show good agreement between 
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theory and the experiments. The Peregrine soliton is the 
first order rational solution of the NTS and the second- 
order rational solution was studied in uni and observed 
experimentally in im. The ratio of maximum ampli¬ 
tude of the rogue wave to the background amplitude is 
3 for the first order rational solution while it is 5 for the 
second order one. Rogue wave solution of the NTS in 
the presence of disorder is also investigated in IfT^ . 

The extension of rogue waves to discrete systems is 
possible. Rogue wave solutions for the discrete NTS, 
Ablowitz-Ladik equation, Salerno model and Hirota 
equation have been presented in ||T31 ITTI [TSl[TH IT71 [TSl 
[m |20l im |22l. We would like to emphasize that the 
ratio of maximum amplitude of the wave to the back¬ 
ground amplitude exceeds 3 and 5 for the first and sec¬ 
ond order discrete solution of the Ablowitz-Ladik lat¬ 
tice, respectively E]. Bludov, Konotop and Akhme¬ 
diev considered the discrete nonlinear Schrodinger 
equation (DNLS) to model an array of nonlinear waveg¬ 
uides f23\ . They constracted a controlled formation of 
a discrete rogue wave by using a proper choice of initial 
field amplitude. In a real experiment, very small degree 
of disorder due to the experimental imperfection always 
exists. In this paper, we show that it has significant role 
on the formation and dynamics of discrete rogue waves 
in an array of nonlinear waveguides. We also predict 
long-living discrete rogue waves. 


2. Discrete Rogue Waves 

The propagation of an optical field in a tight binding 
waveguide array can be described by the following dis- 
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Crete nonlinear Schrodinger equation 

= -JjQVj,, + 'F,_i) + Vj^j + (1) 

where 'Vj is the complex field amplitude at the y'-th 
waveguide, z is the propagation direction, Jj is the cou¬ 
pling coefficient between y'-th waveguide and adjacent 
waveguides, Vj is the propagation constant of the j- 
th waveguide, g is the nonlinear interaction constant. 
We consider attractive interaction, g < 0. As noted in 
El, the solution becomes {-ly'i’je if we replace 
g —> -g. There are two conserved quantities in the sys¬ 
tem. These are total power, P - ^ I'P/I^, and the total 

J 

energy. Our aim is to find discrete rogue wave solu¬ 
tions of the DHLS. Note that a discrete rogue wave is 
the strong localization of the total energy over a few lat¬ 
tice sites II 23 I . 

Suppose the waveguide elements are all identical and 
waveguide separations are the same in the whole sys¬ 
tem. Therefore propagation constant Vj and the cou¬ 
pling coefficient Jj become site independent. In a real 
experiment, the distance between waveguide centers are 
not exactly the same and the waveguide elements are 
not perfectly identical and very small disorder is there¬ 
fore always unavoidable. In this case, the deviations 
of coupling coefficient and propagation constant from 
their site-independent values are very small. We would 
like to emphasize that a small degree of disorder as a 
result of experimental imperfection is always neglected 
on the theoretical study of waveguides. This is reason¬ 
able since almost no contribution comes from such un¬ 
avoidable disorder on fundamental quantum mechani¬ 
cal effects such as Bloch oscillations, Zener tunneling 
and Anderson localization. However we will show be¬ 
low that even very weak disorder changes the physics of 
discrete rogue waves drastically. 

We consider two different types of disorder: diagonal 
and off-diagonal disorders. Off-diagonal disorder is in¬ 
troduced by randomly changing separation between the 
neighboring sites while the sizes of the waveguides are 
all identical. As a result of off-diagonal disorder, the 
coupling coefficient becomes random. Here we assume 
that the coupling constant without disorder equals to 
one, Jo - I- Disorder is introduced by randomly vary¬ 
ing the coupling constant such that Jj e Jo + ei W, where 
ei « 1 is a small dimensionless constant and W are 
random numbers with zero-mean distribution in the in¬ 
terval [-1,1]. The second type of disorder, i.e. diagonal 
disorder, arises due to the random choice of the prop¬ 
agation constant. In a real experiment, there always 
exists very small diagonal disorder due to experimen¬ 


tal imperfections. Therefore, the propagation constant 
Vj varies from one waveguide to another according to 
Vj = Vo H- 62 VT, where To = 0 is set to zero for sim¬ 
plicity, 62 << 1 is a small constant. Since disorder is 
not deliberately introduced, the constants are assumed 
to be very small, 61 ~ 62 < 10^^. We remark that both 
diagonal and off-diagonal disorder are so weak to ob¬ 
serve Anderson localization but they play an important 
role on the formation of the discrete rogue waves. To 
this end, we note that small next neighboring coupling 
is always present in a typical experiment. However, we 
check it has no surprising effect on the formation of dis¬ 
crete rogue waves. So, we neglect higher order cou¬ 
plings in our equation Q. Below we perform numer¬ 
ical computation to study the effect of disorder on the 
formation and dynamics of discrete rogue waves. Be¬ 
fore going further, we define a parameter k as the ratio 
of absolute maximum amplitude the wave can reach to 
the absolute background amplitude. 

Consider first uniform background density. The 
system is known as exhibiting modulational instabil¬ 
ity, which can be investigated by the linearization of 
the dynamical equations around the uniform solution, 
'I'y = where q is the wave number, A is the 

field amplitude and the nonlinear dispersion relation 
reads - 2J cos{q) - gA^. To examine the stability 
of this uniform solution, we add a small perturbation 
by substituting T'j —> 'Vj + d'T j where Q is the 

wave number of the modulation and the dispersion re¬ 
lation can be found by the linearization of DNLS, 
Ql = 87o sin^iQI^) cos(q) {iJq sin^(g/ 2 ) cos(q') + gA^) 
||24l|2l. We conclude that the plane wave solution is 
modulationally unstable whenever becomes nega¬ 
tive. This implies that the unstaggered solution, q - 
0 , is unstable if the nonlinear interaction is focusing, 
g < 0, while the staggered solution, q = tt, is unsta¬ 
ble if the nonlinear interaction is defocusing, g > 0. 
We emphasize that the instability is further enhanced 
in the presence of diagonal and off-diagonal disorder. 
DNLS equation is non-integrable and the above analysis 
is therefore approximate. Furthermore, the linearization 
method shows us the presence of exponentially grow¬ 
ing modes while it gives no predictions about the sub¬ 
sequent stages of evolution. To gain more insight, we 
should perform numerical solution. In the numerical 
simulations, the system is assumed to be subject to pe¬ 
riodic boundary conditions. 

In 12^ . it was discussed that initial conditions must be 
chosen properly to observe discrete rogue waves. We 
first start with a slightly modulated uniform field intro- 
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Figure 1: The absolute of the field amplitude for N = 200, L = 100, 
g = -1 and A = 0.2. The initial form of the field is given by (2). The 
disorder-free case leads to strong localization at the waveguide output 
as shown in (a). The weak off-diagonal disorder with ei = 10 “^ (b) 
and diagonal disorder with ei = 10“^ (c) due to the experimental 
imperfection change the spatial evolution drastically. 


where L is the length of each waveguide in the array and 
A << 1 is the background amplitude. The Fig-1.a plots 
the spatial evolution of the above initial wave. The uni¬ 
form distribution along the lattice is slightly modulated 
at z = 0 and then the modulation grows and reaches its 
maximum at the waveguide output. The total energy is 
strongly localized into a few waveguides at the waveg¬ 
uide output such that the absolute of the maximum am¬ 
plitude is 3 times larger than the background amplitude, 
K - 3. Therefore the discrete wave can be interpreted 
as a discrete rogue wave. Let us now investigate the 
effect of unavoidable weak disorder. Firstly, consider 
weak off-diagonal disorder with ei = 10^^. The pres¬ 
ence of weak off-diagonal disorder leads considerable 
additional energy redistribution among the lattice sites. 
We see that spatial evolution is drastically changed by 
small off-diagonal disorder as can be seen from the Fig- 
1 .b. As opposed to disorder-free case, several peaks oc¬ 
cur at random locations. Note that the ratio k increases 
in one set of random numbers Wj while it decreases in 
another set of random numbers Wj. To conclude, we say 
that if we start with the initial field (|^, what we observe 
in a typical experiment is not like the one in the Fig- 
l.a but in the Fig-l.b. We numerically check that if we 
start with the uniform initial wave 'Tj(z = 0) = Ae^'^ ^ 
instead of we would observe very similar spatial 
evolution as in the Fig-1 .b. It is generally believed that 
a proper choice of initial wave, which is generally chal¬ 
lenging experimentally, is necessary to observe rogue 
waves. Flere we show that even very weak disorder due 
to the experimental imperfection plays more dominant 
role than specific choice of the initial wave on the for¬ 
mation of discrete rogue waves. Let us now consider 
diagonal disorder. The Fig-l.c plots the spatial evolu¬ 
tion of the initial wave 0 in the presence of diagonal 
disorder with 62 = 10^^. As can be seen, very small 
degree of diagonal disorder changes the spatial evolu¬ 
tion similar to the case of off-diagonal disorder. To this 
end, we say that even tiny amount of disorder with ei¬ 
ther ei = 2 X 10^^ or 62 = 2 X 10^^ leads to noticeable 
effects on the spatial evolution. 

Let us now study another type of discrete rogue waves 
and analyze the effect of weak disorder. Consider a pe¬ 
riodical modulation of uniform field parameter 
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where A is the background amplitude, e is a small pa¬ 
rameter and d is the period of the small initial excita¬ 
tion. A question arises. Does modulational instability 
lead to doubly periodic discrete breather ES[ IZTlI (dis¬ 
crete rogue waves that is periodic both along the lattice 
and spatial dimension)? It is not a priori simple problem 
since the spatial evolution sensitively depends on the pa¬ 
rameters. We numerically find that such waves may be 
formed if \g\A^ is not too big. We see doubly periodical 
discrete waves with k - 2.3 and a- = 3 for the parame¬ 
ters A - 0.4 and A — 0.5 when e - 0.2 and g - -0.5. 
If we start with the same parameters but g = -1, dou¬ 
bly periodic structure would be lost. The Fig-2.a plots 
doubly periodic discrete wave when A - 0.4. The evo¬ 
lution starts with a periodically modulated background 
and the modulation increases until it reaches its max¬ 
imum. Then the system returns to the original back¬ 
ground. This repeats periodically in spatial dimension. 
Let us now study the effect of weak disorder due to 
the experimental imperfections. The Fig-2.b and Fig- 
2.C plot spatial evolution in the case of off-diagonal and 
diagonal disorders, respectively. The modulational in¬ 
stability as a result of small degree of disorder is slower 
and has almost no effect when z < 50. Therefore the first 
cycle remains almost unaffected as shown in the figures 
2.b and 2.c. The disorder-related modulational insta¬ 
bility dominates at later times and the spatial evolution 
becomes chaotic. As a result, the wave is not doubly pe¬ 
riodic discrete wave anymore but a discrete rogue wave. 
The presence of disorder increases k from 2.3 to 3 with 
our set of random numbers Wj. 

We have shown that weak disorder plays a vital role 
on the formation and spatial evolution of discrete rogue 
waves. Furthermore, an experimental realization of ini¬ 
tial waves like is generally challenging. Therefore, 
we propose to start with the uniform initial density and 
let the weak disorder due to the experimental imperfec¬ 
tions lead to discrete rogue waves in the system. One 
can then experimentally observe spatial evolution like 
the one as in the Fig-l.b and Fig-l.c. Now, a ques¬ 
tion arises. What is the long time behavior of the dis¬ 
crete rogue waves? One intuitively expects that a cou¬ 
ple of peaks appears at random lattice sites and disap¬ 
pears shortly and then another peaks appear in other lat¬ 
tice sites and disappear shortly too. In the same time, 
chaotic fluctuations happen on the background. As a re¬ 
sult, one expects that the positions of the peaks change 
(not continuously) along the propagation direction. Let 
us now answer the question. If z is not very large, the 
dynamics is as what it is expected. Discrete rogue waves 
with K ~ 3 appear at random positions and they decay 
due to the coupling to the neighboring sites. If z is large. 


there are more peaks in the system and some of them 
are close enough to each other to start collision. Some 
collisions in the nonlinear medium excite discrete rogue 
waves with k ~ 5. The numerical calculations reveal 
that the discrete rogue waves with k ~ 5 don’t decay 
to the adjacent sites. In this way, long-living discrete 
rogue waves are formed. The Fig-3.a displays long dis¬ 
tance spatial evolution of the initially uniform field. As 
can be seen from the figure, the amplitudes of the three 
peaks change slightly with z while dynamical changes 
occur chaotically at other lattice sites. The peaks appear 
from nowhere and almost maintain their profiles. The 
second message of this paper is the prediction of long- 
living discrete rogue waves. Let us now understand the 
origin of such waves. The standard perturbation method 
fails here since such waves with k ~ 5 appear at large z 
and the system is also chaotic. To understand the for¬ 
mation of long-living rogue wave, consider an initial 
peak at the middle of the lattice with N = 200. There¬ 
fore we suppose the initial field reads '¥j(z = 0) = A for 
j + 100 and T'j(z = 0) = pA for j - Njl - 100, where 
p > 1 is a constant and A = 0.4. We expect that the ex¬ 
cess energy at j - 100 starts to diffract to adjacent sites 
unless nonlinear interaction is strong enough (if gp^A^ 
is not larger than Jq). If p is large enough (if p > 5 for 
our numerical parameters) then nonlinear interaction is 
dominant and diffraction affects perturbatively. In this 
case, one can readily say that the amplitude of the peak 
stays almost the same in the early stage of dynamics. At 
around z = 50, a couple of peaks start to appear in the 
system due to modulational instability. This increases 
the effect of diffraction on adjacent sites of the peaks. 
Another question is now arising. Is the peak at j - 100 
destroyed at a large value of z? The answer is no as one 
can see in the Fig.3.b. This can be understood as fol¬ 
lows. The coupling of the site j = 100 to the adjacent 
sites can be enhanced by the existence of another peak 
at the adjacent sites, which can be formed by modula¬ 
tional instability. It is well known that rogue wave is a 
rare event and the probability of rogue wave formation 
exactly at the adjacent sites is very small. Therefore the 
peak at j - 100 is long-living and can be called a lo¬ 
calized discrete soliton on a finite background. Having 
studied the dynamics of the single initial peak, let us 
now analyze our original problem. Our original system 
is initially uniform and small disorder with 62 — 10“^ 
leads to modulational instability. Rogue waves with 
K ~ 3 appear at around z = 100 as we have already seen 
in the Fig. 1 .c. There are more than one peak in the sys¬ 
tem and collisions among them lead to higher amplitude 
peaks. At around z ~ 300, some peaks with k ~ 5 ap¬ 
pear. As discussed above, they survive very long since 
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their couplings to the adjacent lattice sites are negligi¬ 
ble. Additionally, no collisions occurs between these 
peaks since they don’t tunnel to adjacent sites. This pre¬ 
vents the formation of higher order rogue waves with 
K » 5. This is the reason why k starts from one and in¬ 
creases up to nearly 5 in our numerical simulation. We 
call it long-living discrete rogue wave since the high- 
amplitude peaks appear from nowhere and don’t disap¬ 
pear. 

3. Conclusion 

In this paper, we have studied discrete rogue waves in 
an array of waveguides. We have discussed that small 
degree of disorder always exists in a typical experiment. 
We have considered diagonal and off-diagonal disor¬ 
ders. There are two main results of this paper. Firstly, 
we find that small degree of disorder leads to significant 
effects on the formation and dynamics of discrete rogue 
waves. This is interesting since such a weak disorder 
has almost no role on Bloch oscillation, Zener tunneling 
and Anderson localization. Secondly, we predict long- 
living discrete rogue waves. The long-living discrete 
rogue wave studied here appears from nowhere and al¬ 
most maintains its profile instead of vanishing without 
trace. 
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Figure 2: The absolute of the wave function for N = 200, g = -0.5, 
A = 0.4, 6 = 0.2 and d = 20. The initial form of the field is given by 
(3). A doubly periodic discrete rogue wave is formed for the disorder- 
free case as shown in (a). The doubly periodic character is lost in 
the presence of weak off-diagonal disorder with e\ = 10 “^ (b) and 
diagonal disorder with 62 = 10 “^. 


Figure 3: The uniform initial field, = 0) = A, leads to long-living 
discrete rogue waves as can been from the Fig. a. We plot the figure 
for 0 < z < 6000 but the long-living waves survives for a much larger 
z. In the figure (b), a single site initial peak on the uniform back¬ 
ground is considered, ^j(z = 0) = A for j 100 and Ty(z = 0) = 5A 
for j = 100. The fig.b show us that if a peak is large enough than the 
background amplitude, then nonlinear interaction is more dominant 
than diffraction. Note that the long-living peaks occur due to modula- 
tional instability at large z- The parameters in the figures are given by 
^2 = 10^^, g = -1, W = 200 and A = 0.4. 
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